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Patterning is a familiar approach for imparting novel functionalities to free surfaces. We extend the 
patterning paradigm to interfaces between crystalline solids. Many interfaces have non-uniform internal 
structures comprised of misfit dislocations, which in turn govern interface properties. We develop and 
validate a computational strategy for designing interfaces with controlled misfit dislocation patterns by 
tailoring interface crystallography and composition. Our approach relies on a novel method for predicting 
the internal structure of interfaces: rather than obtaining it from resource-intensive atomistic simulations, 
we compute it using an efficient reduced order model based on anisotropic elasticity theory. Moreover, our 
strategy incorporates interface synthesis as a constraint on the design process. As an illustration, we apply 
our approach to the design of interfaces with rapid, 1-D point defect diffusion. Patterned interfaces may be 
integrated into the microstructure of composite materials, markedly improving performance. 

Composite materials with exceptional properties, such as radiation resistance\ extreme thermo -mechanical 
stability^, or high strength and ductility^ may be created by controlling the interfaces between composite 
constituents, in addition to the constituents themselves. One way to customize interface properties is by 
tailoring their internal structure. However, unlike free surfaces, which may be directly patterned'*"^, interface 
structure is difficult to control because interfaces are buried within the material and cannot be directly accessed. 

We propose an approach for computational design of patterned interfaces by tuning their crystallography and 
composition. We focus specifically on interfaces of the type illustrated in Fig. 1, whose internal structure contains 
coherent regions separated by inherent line defects, known as misfit dislocations^^"^^. Interface crystallography 
and composition are the natural "design space" for such interfaces because they govern the misfit dislocation 
pattern. This pattern, in turn, determines interface properties such as impurity precipitation^^"^^, point defect 
mobility^^"^^ and shearing resistance^^'^°. Thus, given the misfit dislocation pattern that yields a desired property, 
our goal is to predict the interface crystallography and composition that generate that pattern and — ideally — 
specify a synthesis route for that interface. 

Misfit dislocation patterns may be predicted using atomic-level interface models^^'^^'^^. We view this approach 
as a well-posed "forward problem." By contrast, interface design— i.e., finding the interface that yields a desired 
misfit dislocation pattern — is an "inverse problem"^^. In principle, it may be solved by repeatedly executing the 
forward problem over the design space. In practice, such a procedure is not feasible because of the high dimen- 
sionality of the design space (discussed below) and the substantial effort required for each forward problem 
solution. 

Our approach to interface design is to construct a mesoscale (as opposed to atomic-level) model that predicts 
misfit dislocation patterns with accuracy comparable to atomistic simulations, but at a fraction of the cost. Our 
model is a reduced order model because it replaces the millions of variables associated with atomic positions with 
<15 variables needed to describe misfit dislocations. This reduced order mesoscale model, or ROMM, enables 
rapid scanning over the interface design space and thereby opens the path to computational design of patterned 
interfaces. To demonstrate its effectiveness, we use our ROMM to design interfaces that exhibit rapid, 1-D 
diffusion of point defects. 

Methods 

Design space. Solid-state interfaces exhibit limitless variety We focus on the subset of atomically sharp interfaces between metals 
whose crystal structures may be mapped into each other by uniform deformation gradients. Far from being restrictive, this subset defines a 
high-dimensional design space that includes all stress-free and uniformly strained face- and body-centered cubic (FCC, BCC) and tetragonal 
(FCT, BCT) metal pairs. Grain boundaries also fall within this subset. Extensions of our ROMM to other types of interfaces will be addressed 
in the Discussion section. 

We describe an interface's crystallography by comparing it to a single crystal reference state^^'^^. Any interface of the type studied here may 
be described by two uniform deformation gradients: aF and bF, applied respectively to the parts of the reference crystal above (A) and below 
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Figure 1 | Ag terminal plane in a Ag/V interface (crystallography as in 
Results section with 9=1°). Interface atoms are colored by their energy, 
showing coherent regions separated by misfit dislocations. 

(B) the X-Z plane, as shown in Fig. 2. This mapping generates the "natural" interface 
state, with the interface along the x-z plane. aF and bF are chosen such that, in the 
natural state, parts A and B match the prescribed crystal structures and orientations of 
the solids joined at the interface. Misfit dislocation patterns, such as those in Fig. 1, are 
generated once the interface in its natural state is relaxed. 

aF and bF may be represented as 3 X 3 matrices^*^'^''. Together, they contain 
eighteen components. Six of them correspond to rotations of A and B. Since equal 
rotations of A and B along the interface plane do not affect interface structure, only 
five rotation components are needed to describe interface crystallography. The 
remaining twelve components of aF and bF describe strains. However, only the 
relative strain between A and B affects interface structure, reducing the number of 
independent strains to six. Therefore, there are in total eleven crystallographic degrees 
of freedom (DOFs) available for interface design. For grain boundaries, the strains are 
fixed, reducing the DOFs to the well-known five^^. 

In addition to crystallographic DOFs, interface structure also depends on the 
elastic properties of the adjoining materials^^'^^. A solid may have between two and 
twenty-one independent elastic constants, depending on its symmetry. Thus, a 
complete description of interface elasticity requires at least four elastic constants (if 
the adjoining materials are isotropic) and as many as forty-two. In practice, elastic 
constants depend on composition and crystal structure. Thus, rather than enumerate 
all their possible combinations, we simply refer to them as a potentially very large set 
of "compositional" DOFs available for interface patterning. 

State variables. Given the crystallographic and compositional DOFs of an interface, 
the resulting misfit dislocation pattern may be determined using atomistic 
simulations. Such simulations are costly because they track the positions of all atoms 
near the interface. It is far more efficient to describe interface structure in terms of the 
misfit dislocations themselves, instead of atoms. Misfit dislocation patterns for the 
interfaces of interest here consist of m sets of misfit dislocations. Each set contains an 
infinite number of straight, parallel, equally spaced dislocations, as illustrated in 
Fig. 3.b)-d). Thus, a minimal description of each set requires the dislocations' line 
direction their spacing di, and Burgers vector hi (i= 1 . . .m). 

is a unit vector in the interface plane and so may be specified using a single 
variable while di and hi require one and three variables, respectively. Therefore, as few 
as 5m state variables are required to describe all m sets of dislocations, as opposed to 
the millions of variables needed to specify locations of all atoms near an interface. 
Potential additional state variables will be considered in the Discussion section. 

Reduced order mesoscale model (ROMM). We construct a ROMM that predicts, 
without recourse to atomistic simulations, the 5m state variables of misfit dislocation 




Figure 2 | Interface crystallographic DOFs may be described with 
deformation gradients aF and bF, which map from a single crystal 
reference state to the natural state. 



a) Reference state b) case 1 




c) case 2 d) case 3 




Figure 3 | (a) Atom positions at the interface plane in the coherent 
reference state for a Cu/Nb interface with 6=1°. (b) "Case 1" is a solution 
of the qFBE (Eqn. 1) that uses Burgers vectors hi and ^2 shown in (a), (b) 
"Case 2" uses hi and h^ while (c) "case 3" uses ^2 and h^. di and denote 
misfit dislocation spacings and directions. 

patterns from the 15+ DOFs describing interface crystallography and composition. 
The crystallographic DOFs and state variables are related by the quantized Frank- 
Bilby equation (qFBE)^^'^''-^^: 

(AF-'-BF-)i5=E"^,(^f)(;, (1) 

The left hand side is derived from a Burgers circuit enclosing a probe vector p that lies 
in the interface plane. It expresses the sum of misfit dislocation Burgers vectors 
crossed byp (averaged over all possible origins of p) from deformation gradients aF 
and bF. The right hand side of the qFBE writes out this sum in terms of up to m = 3 sets 
of misfit dislocations. 

Solutions for ^i and di maybe obtained from the qFBE using O-lattice theory^^'^^'^^. 
The determination of hi involves a subtle problem in anisotropic elasticity^*, which we 
recently solved^^. However, these two steps do not yield a unique solution. Rather, 
they provide as many candidate solutions as there are combinations of linearly 
independent Burgers vectors in the reference state^''-^^ We determine which of these 
candidate solutions is the correct one by calculating their interface energies, y> and 
choosing the one with lowest y. 

Misfit dislocations may be viewed as Volterra dislocations that have been inserted 
into the coherent reference state^^, suggesting that y be expressed as 

y=ycoh+yeiastic+y core + y relax ~r • • • (2) 

Ycoh is the interface energy of the coherent reference state and the remaining terms are 
due to misfit dislocations: elastic strain energy Ydastio core energy Ycore> energy due to 
relaxations of the misfit dislocation network Yreiax> and perhaps additional terms that 
have not been identified. 

For our purposes, it is not necessary to calculate the absolute value of y, but rather 
only differences in y between the candidate solutions of the qFBE. Ycoh does not 
contribute to these differences because each solution represents different Volterra 
dislocations inserted into the same coherent reference state. We further make the 
hypothesis— to be validated later— that core and relaxation energies are not the main 
contributors to differences in y. We therefore concentrate on computing the elastic 
energy contribution: Yeiastic- 

After enumerating all the solutions to the qFBE (Eqn. 1), we compute the complete 
elastic fields for each one, as described in Ref 27. These fields satisfy four boundary 
conditions: (i) absence of net far-field strains, (ii) consistency of far-field rotations 
with the misorientation prescribed in aF and bF, (iii) no net tractions along the 
interface, and (iv) a position-dependent displacement discontinuity, i.e. "disregis- 
try"^^, across the interface consistent with the misfit dislocation pattern of interest. 
Disregistries may be compared directly with atomistic simulations^'''^"'^^'^*^. Our cal- 
culation is fully anisotropic and may be performed for any number of independent 
elastic constants. It is therefore capable of accounting for all compositional DOFs 
under consideration here. 
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Figure 4 | Interface energies computed as a function of 8 using our ROMM and atomistic modeling for (a) Cu/Nb and (c) Ag/V. Filled circles indicate 
atomic models whose disregistry was analyzed. The ringed numbers next to them state the case that best matches the atomic disregistry. Ad hoc 
parameters P, Q, and R (Eqn. 4, 5, and 6) for (b) Cu/Nb and d) Ag/V. 



When all boundary conditions are satisfied, the interface is free of far-field stresses. 
However, short-range stresses persist near the interface. These short-range fields are 
used to compute interface elastic energies: 

yelastic=^ // (Tshortn.'AUdis dS. (3) 
JJ A(ro) 

A is the domain of integration: a single unit cell of the misfit dislocation pattern. To 
eliminate the unphysical divergence of stresses at dislocation cores, A excludes the 

region within a cutoff distance ro = l/4|5| of the dislocation cores. Oshonn and Aw dis 

are the traction and disregistry vectors at the interface, respectively. By applying Eqn. 
3 to all candidate solutions of the qFBE (Eqn. 1), we predict which one has lowest 
Yeiastic and therefore describes the likeliest misfit dislocation pattern. 

Results 

Validation. Our ROMM is a significant simplification of a complex 
problem and requires validation. To this end, we compare the output 
of our ROMM with atomistic calculations. We focus on heterophase 
interfaces between {111} planes of FCC metals and {110} planes of 
BCC metals. Such interfaces contain misfit dislocations of unlike 
character and asymmetric arrangement^^'^^'^^'^^. Thus, they provide 
an opportunity for rigorous validation of our ROMM. They are also 
convenient for atomistic simulations because embedded atom 
method (EAM)^^ potentials are available for several FCC/BCC 
binaries. 

We validate our ROMM against four interface compositions: Cu/ 
Nb^°, Ag/V'*\ Cu/Fe^^ and Cu/Mo^l These choices fix the elastic 
constants, crystal structures, and lattice parameters of the adjoining 
constituents. Because we restrict attention to interfaces along FCC 
{111} and BCC {110} planes, only one crystallographic DOF remains 



to be specified: the twist angle 0 describing the relative rotation of the 
crystals parallel to the interface plane. 

We measure 0 with respect to the Nishiyama- Wasserman orienta- 
tion relation, where a BCC <100> direction is parallel to a FCC 
<110> direction (see supplementary note)^^'^^ 0 = n/?> — 
cos~^ (l / Vs) ~ 5.26° yields the Kurdjumov-Sachs orientation rela- 
tion (KS OR)^^. Due to the symmetry of the interface planes, all 
crystallographically distinct interfaces fall within 0° < 0 < 15°. 
However, we limit our analysis toO°<0<lO° because for greater 
twists, misfit dislocations are too closely spaced to characterize reli- 
ably in atomic models. 

For any composition and 0, the qFBE has three distinct candidate 
solutions^^. Each consists of two sets of misfit dislocations (m=2) 
and corresponds to one of three combinations of interfacial Burgers 
vectors, defined in the coherent reference state. For example. Fig. 3. a) 
shows the Burgers vectors in the reference state of a Cu/Nb interface 
with 0=1°. The first candidate, termed "case 1," uses Burgers vectors 
hi and ^2- "Case 2" and "case 3" use Burgers vectors hi, h^ and h2, ^3, 
respectively. Each case yields a different misfit dislocation pattern, as 
shown in Fig. 3.b) -d). The coherent reference state is the same for all 
three cases. aF and bF are chosen such that hi, h2, and ^3 map to 

-(110) -type FCC directions, hi maps to a <100>-type BCC dir- 
ection, and h2, and I73 map to l/2(lll)-type BCC directions. 

Using our ROMM, we compute Yeiastic of all three cases for each 
composition and 0 of interest. For all interfaces, we also construct 
atomic- scale models by joining cylindrical FCC and BCC blocks 
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following the required interface crystallography. The models are 
large enough to contain a representative area of the misfit dislocation 
pattern and to avoid elastic images from free surfaces. The models are 
relaxed using EAM potentials^^"^^ and analyzed to extract interface 
energies and misfit dislocation structures (see supplementary note). 
All ROMM calculations use lattice parameters and elastic constants 
predicted by these EAM potentials. 

Fig. 4.a) compares Yekstic i^om our ROMM with y from atomistic 
simulations for Cu/Nb interfaces. Since we are only interested in the 
relative energies of the three cases, both the ROMM and atomistic 
data are shifted so that their energy minima occur at 0 J/m^. The 
ROMM predicts that case 3 has lowest Yekstic for Furthermore, 
Yekstic for case 3 is in near perfect quantitative agreement with y for 
0>1°. 

Fig. 4.c) shows a similar comparison for Ag/V interfaces. Here, the 
ROMM predicts that case 1 has lowest Yekstic for all 0 outside 
4.25°<0<5.25°, where 

Yekstic is lowest for case 2. Yekstic arid Y are 
in qualitative agreement over the entire twist angle range and in 
quantitative agreement for 0> 5°. As detailed in the supplementary 
note, we find comparable agreement between the ROMM and ato- 
mistic interface energies for the remaining two compositions we 
studied. 

Agreement between Yekstic and y is not sufficient to validate our 
ROMM. For that, we must determine whether the lowest energy 
cases predicted by the ROMM match the misfit dislocation patterns 
in our atomistic simulations. To this end, we compare the disregis- 
tries of the two models (see supplementary note). Each of the three 
qFBE solutions predicts a different misfit dislocation pattern and 
therefore also a different disregistry. Our goal is to compare the 
disregistries of all three cases with that found in atomistic simula- 
tions. Our ROMM is validated if the case with lowest Yekstic has the 
best match with the atomistic disregistry. 

As shown in Fig. 4.a) and c), the disregistry analysis is in agree- 
ment with our ROMM predictions for all Cu/Nb and Ag/V interfaces 
except Cu/Nb at 0=0°. We attribute the disagreement to the recon- 
struction of the misfit dislocation network that is known to occur at 
that interface^^. One further case of disagreement where dislocation 
network reconstruction occurs is found for Cu/Mo at 0 = 0° (see 
supplementary note). However, the agreement between our 
ROMM and the atomistic models is excellent, overall. 

Our ROMM may be compared with several ad hoc parameters 
proposed previously to determine which of the cases predicted by 
the qFBE is likeliest. BoUmann suggested that the likeliest case mini- 
mizes^^ 



(4) 



The rationale for this parameter is analogous to the "Frank rule" for 
predicting dislocation reactions^^. Similarly, Ecob and Ralph propose 
to distinguish between cases using parameters'^ 



didj 



or 



(5) 



(6) 



Fig. 4.b) and d) plot these parameters for Cu/Nb and Ag/V interfaces. 
Comparing with Fig. 4.a) and c), we see that none of them predicts 
the misfit dislocation patterns seen in atomistic models. For example, 
for Cu/Nb, all three parameters favor case 2, while the true interface 
structure is case 3. We therefore view our ROMM as superior to these 
parameters and as validated for the purpose of computational design 
of patterned interfaces. 



Application to interface patterning. A validated ROMM for pre- 
dicting misfit dislocation patterns may be used for computational 
interface design. Potential applications include interface patterning 
for templated precipitation of impurities^^'^^'^^, tailored shearing 
resistances^, or twin nucleation^^. As a demonstration, we apply 
our ROMM to the design of interfaces with fast, 1-D point defect 
diffusion. 

Previously, we detailed a distinctive mechanism of point defect 
diffusion at a Cu/Nb interface^^"^^. In it, interfacial vacancies or inter- 
stitials undergo rapid migration parallel to one set of misfit disloca- 
tions and considerably slower migration perpendicular to it. This 
mechanism may also operate at other interfaces provided that their 
misfit dislocation pattern meets two criteria. First, it must contain 
one set of misfit dislocations of nearly perfect screw character and 
another of non- screw character. Second, the distance between misfit 
dislocation intersections along the fast migration direction must be 
small, i.e. below ~2 nm. 

We use our ROMM to design FCC/BCC interfaces that meet these 
two criteria and therefore exhibit rapid 1-D point defect diffusion. 
We consider interfaces of the type described in Validation subsec- 
tion, but rather than restricting ourselves to pure elements, we use 
alloys of Cu and Ni for the FCC side and alloys of Nb and V for the 
BCC side. By changing the alloy composition, we may continuously 
vary the lattice parameters and elastic constants of the adjacent crys- 
tals. We assume that all of these properties follow Vegard's law, i.e. 
that they may be interpolated linearly in composition between the 
single -component properties. 

For example, the lattice parameter of the Cu-Ni alloy may be 
expressed as 



^CuNi — ^Cu ( 1 — ^m) + i^Ni^h 



(7) 



where d^i is the atomic fraction of Ni in the alloy. Thus, d^i is the only 
compositional DOF for the FCC component. Similarly, the atomic 
fraction of V, Sy, is the only compositional DOF for the BCC com- 
ponent. Replacing Vegard's law with a nonlinear composition-prop- 
erty relation would not increase the number of compositional DOFs 
and so would not change the interface design process, although it 
might alter its outcome. The only crystallographic DOF we will con- 
sider is the twist angle 0. 

Using our ROMM, we find that the design space con- 

tains regions with case 2 and case 3 misfit dislocation patters, as 
shown in Fig. 5. a). With this information, we plot envelopes where 
the criteria for rapid 1-D diffusion are met. Inside the red, airfoil- 
shaped regions in Fig. 5. a), one set of dislocations is within 5° of 
perfect screw character while the other set of dislocations is non- 
screw. Within the region shaded gray, the spacing between misfit 
dislocation intersections along the screw set of dislocations is smaller 
than 2 nm. The intersection between the red airfoils and the gray 
shaded region represents misfit dislocation patterns that meet the 
criteria for rapid, 1-D diffusion. 

Ideally, a design for a patterned interface identifies not only crys- 
tallographic and compositional DOFs, but also provides a synthesis 
route for that interface. In this example, we may incorporate syn- 
thesis as an additional constraint on the design process. For example, 
interfaces between FCC {111} and BCC {110} planes may be made 
using physical vapor deposition (PVD)^^. PVD gives interfaces in the 
KS OR, where 0—5.26°. Thus, the design envelope with a PVD syn- 
thesis route is obtained by restricting the 3-D design space to a planar 
submanifold with 0—5.26°, shown in Fig. 5.b). Other interface pro- 
cessing methods, such as molecular beam epitaxy (MBE)'*^, diffusion 
bonding^°, precipitation^\ or severe plastic deformation (SPD)^^, 
may be represented by different constraints on the design space. 

In principle, the interface design process described above may be 
carried out using atomistic modeling. In practice, the resources 
required would be prohibitive. Nearly 69 thousand distinct interfaces 
were examined to create the 3-D map in Fig. 5. a). Assuming an 
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Figure 5 | (a) The 3-D 5Ni-5v-6 design space and (b) a planar section through it for 9—5.26°. Domains corresponding to case 2 and case 3 are deUmited by 
a 2-D surface in (a). The red, airfoil- shaped envelopes labeled "(|) < 5°" enclose interfaces that have one set of nearly screw misfit dislocations as well as one 
non-screw set. The shaded region labeled "s < 2 nm" encloses interfaces with misfit dislocation intersections less than 2 nm apart along the set of 
misfit dislocations with screw character. The solution space for interfaces with rapid 1-D point defect diffusion is the intersection between this region and 
the red airfoils. 



interatomic potential is available, we estimate 5 person-hours and 24 
CPU hours to compute one interface structure using atomistic mod- 
eling. If a potential is not available, then considerable additional 
human effort is required to create one. Alternatively, much greater 
computing resources may be used to carry out first principles ato- 
mistic simulations. Thus, constructing the map in Fig. 5. a) using 
atomistic modeling is infeasible. By contrast, <7 CPU minutes are 
required to compute an interface structure using our ROMM. 
Furthermore, since the ROMM takes DOFs such as elastic constants 
as direct inputs, it may be applied to a wide range of compositions 
without further effort. 



Discussion 

The promise of using designer solid-state interfaces to impart unpre- 
cedented new functionalities to materials has been widely recog- 
nized^"^'^^. We presented a strategy for rapid computational design 
of interfaces with tailored misfit dislocation patterns. Computational 
design of patterned interfaces may be viewed as an inverse problem 
where the DOFs that yield an interface with prescribed state variables 
are found. Our strategy for solving such inverse problems relies on a 
reduced order mesoscale model, or ROMM, for predicting interface 
state variables from interface DOFs. We validated our ROMM 
through direct comparisons with atomistic models and applied it 
to an illustrative interface design problem. Synthesis routes may be 
incorporated as constraints in the interface design process. 

We illustrated our ROMM-based strategy on the example of inter- 
face design for fast, 1-D point defect diffusion. A critical piece of 
knowledge that enabled that design was the quantitative description 
of the mechanism of point defect diffusion at Cu/Nb interfaces, 
which was already available^^"^^. Such a model linking misfit disloca- 
tion patterns to interface properties may itself be viewed as a kind of 
ROMM. Further applications of interface design will also require 
quantitative models linking interface structure to interface prop- 
erties. Several are already available. For example, interfacial precipi- 
tates wet high-energy regions near misfit dislocations while avoiding 
low-energy, non-wetting coherent patches^^. Thus, control of inter- 
face structure may be used for templated precipitation^^, in analogy 
to free surfaces patterned to control water nucleation^ or to enable 
self-cleaning^. Similarly, impurity trapping^^'^'* and shear resistance^° 
differ for coherent patches and misfit dislocations, opening a path to 
controlling interfacial segregation and mechanical properties. 



Although motivated by computational design of patterned inter- 
faces, the ROMM described here may also find applications in pre- 
dicting phase transformations or microstructure evolution. These 
processes depend on interface structure and energy^^'^^, both of 
which may be computed using our ROMM. For example, Cu/Nb 
interface energy is minimized at 0—2° in Fig. 4.a). Thus, given a 
Cu/Nb interface initially at some other 6 and assuming no con- 
straints from surrounding grains, we expect the adjoining crystals 
to rotate until the minimum energy misorientation is reached. This 
reasoning may be extended to predict orientations of martensite 
habit planes^^ and grain boundary energies in polycrystals^^. 

Our ROMM distinguishes between candidate interface structures 
predicted by the qFBE based on strain energies Yeiastk calculated 
using anisotropic elasticity theory^^'^^. This is the underlying reason 
for the better performance of our ROMM compared to previously 
proposed parameters^^'^^ It is nevertheless remarkable that Yeiastic is 
often in quantitative agreement with y computed using atomistics. 
The true interface energy, y, involves several terms besides yeiastic^ ^-g- 
contributions from dislocation network relaxations and from dis- 
location core energies. Indeed, the former appear to play a role when 
the agreement between yeiastic ^nd y is not perfect. For example, the 
misfit dislocation network in Cu/Nb at 0 = 0° is known to reconstruct 
into a hexagonal pattern^^, explaining why y<yeiastic in Fig. 4.a). We 
have found evidence for similar reconstructions in Ag/V for 0<4.5°. 

More surprisingly, the agreement between yeiastic and y is very 
good despite the low spacing between misfit dislocations in the 
FCC/BCC interfaces in the Results section: in some cases, it is com- 
parable to the dislocation core size, i.e. < 1 nm. The fact that yeiastic 
nevertheless provides a good estimate for y is likely not because 
dislocation core contributions ycore ^1^7 be neglected, but rather 
because they are related to the elastic energy contribution. Such a 
relationship may be rationalized through the Peierls-Nabarro 
modeP^'^° or analogous arguments^^'^^. 

Numerous extensions of the ROMM described here may be envi- 
sioned, for example to interfaces with three sets of misfit dislocations, 
interfaces whose reference states are not single crystals (e.g. grain 
boundaries with small misorientations from coherent states), or 
interfaces between crystals of finite size (to account for free surfaces 
or layered morphologies). By increasing the number of state vari- 
ables, the ROMM may be extended to account for misfit dislocation 
network relaxations, dislocation core structures, solute atmospheres, 
or non-equilibrium interface states. 
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Expanding the subset of interfaces treated by the ROMM would 
require increasing the number of DOFs it takes as inputs. For 
example, modeling interfaces between crystals that cannot be 
mapped to each other using uniform deformation gradients alone 
may require new DOFs for sublattice atom shuffles. Additional DOFs 
may be required to describe non-elastic energy contributions, e.g. 
electrostatic ones in interfaces between ionic solids. Line defects 
other than dislocations may be important for some interfaces^^'^^. 

Since the time of Read and Shockley^'*, it has been argued that 
misfit dislocation models of interfaces must break down when the 
distance between dislocations becomes small. However, in view of 
the numerous possible improvements and extensions mentioned 
above, it is not clear at what dislocation spacing this breakdown 
would occur or how far it may be pushed back. As ROMMs for 
interface structure become more complex, determining the limits 
of their validity may increasingly require recourse to formal uncer- 
tainty quantification (UQ) methods, such as ones based on Bayesian 
inference^^. 
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